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ABSTRACT 

The synchrotron reflection scenario recently proposed to explain 7-ray flares 
observed from blazar jets is studied. Our analysis takes into account the angular 
distribution of the beamed radiation, the finite extent of the scattering region, and 
light travel-time effects. We compare energy densities and powers for synchrotron, 
SSC, reflected synchrotron (RSy), and external Compton (EC) scattering processes. If 
the width of the scattering layer is much larger than TR' B , where T and R' B denote the 
bulk Lorentz factor and comoving-frame radius of the plasma blob, respectively, then 
the ratio of the RSy and synchrotron energy densities ~ 4 r 3 riBLR cr T^B > where tiblr 
is the mean particle density in the broad line region (BLR). Our results imply that 
Thomson-thick scattering regions of narrow extent must be present for the synchrotron 
reflection mechanism to operate effectively. This process seems unlikely to cause flares 
in lineless BL Lac sources, where X-ray and TeV flares are common and the BLR is 
thought to be weak or absent. We sketch time profiles of flares for various scenarios, 
including a model where the blob is energized by sweeping up surrounding material. 

Subject headings: galaxies: active — galaxies: jets — radiation mechanisms: 
nonthermal 

1. Introduction 

More than 50 blazar-type AGNs have been detected with high confidence by EGRET to emit 
7-rays above 100 MeV (Matox et al. 1997). These sources are identified with flat-spectrum radio 
sources classified as BL-Lac objects or quasars. Many of these objects exhibit variability on all 
wavelengths, with some of the most rapid variability, on time scales of hours to days, observed at 
the highest 7-ray energies (e.g., Bloom et al. 1997; Wagner et al. 1995; Mukherjee et al. 1997). 

The large apparent luminosities in combination with the short variability time scales provide 
evidence for the widely accepted relativistic jet model for AGNs (for recent reviews, see Schlickeiser 
1996 and Hartman et al. 1997), according to which the radio-7-ray emission from blazars is 
emitted via nonthermal synchrotron radiation and Comptonization of soft photons by energetic 
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particles in relativistic outflows powered by accreting supermassive black holes. Soft photons 
which are Compton-scattered to produce the 7-ray emission include internal synchrotron photons 
(e.g., Marscher & Gear 1985, Maraschi et al. 1992, Bloom & Marscher 1996) and accretion-disk 
radiation which enters the jet directly (Dermer & Schlickeiser 1993) and after being scattered by 
surrounding BLR clouds and circumnuclear debris (e.g., Sikora, Begelman Sz Rees 1994; Blandford 
& Levinson 1995; Dermer, Sturner, & Schlickeiser 1997; Protheroe & Biermann 1997). 

It has recently been proposed (Ghisellini & Madau 1996; hereafter GM96) that the beamed 
synchrotron radiation, after scattering off a cloud near the jet trajectory and reentering the jet, 
can be a source of copious soft photons and lead to a pronounced flare of very short duration 
as the relativistic jet plasma passes through the cloud. Wehrle et al. (1998) argue that this 
mechanism might explain the February 1996 7-ray flare observed from 3C 279. A detailed analysis 
which correctly accounts for causality effects and the finite width of the scattering layer was not 
performed by GM96, and such a treatement is required before blazar flare spectra and light curves 
can be modeled. 

Here we examine this model in more detail for a simple geometry of the BLR in the limiting 
regime where the nonthermal jet electron distribution does not evolve. In §2 we describe the 
model. Numerical calculations of the magnetic-field and photon energy densities and synchrotron 
and Compton powers are presented in §3. Application to blazar flares is made and the effects of 
blob energization on time profiles of flares are indicated. We summarize in §4. 

2. Description of the Model 

The idealized system, illustrated in Figure 1, is assumed to have the following qualities: A 
plasma blob, which in the comoving frame is spherical with radius R' B ~ 10 15 - 10 cm, moves 
with bulk Lorentz factor T ~ 10 and velocity (3rc along the z axis which defines the jet. The 
blob is filled with relativistic, nonthermal electrons and pairs which are isotropically distributed 
in the comoving frame with a comoving power-law Lorentz factor distribution n e {^) oc 7 _p for 
7i < 7 < 72- The nonthermal lepton density in the jet is denoted by n e j e tj and a randomly 
oriented magnetic field B, in equipartition with the energy in relativistic leptons, is present in 
the blob. An accretion disk is located at z = and radiates with a luminosity ~ 10 - 
10 46 erg s _1 . The disk spectrum is approximated by a thermal blackbody of ~ 10 5 K peaking 
at the characteristic maximum energy of a standard Shakura-Sunyaev disk spectrum (Shakura & 
Sunyaev 1973) for a central black hole of mass M ~ 10 8 M Q . 

The accretion-disk/jet system is assumed to be surrounded by a spherically distributed 
layer of scattering material which represents the BLR. This layer extends from r- m to r out and 
has a radial Thomson depth tt,blr ~ 0.01 - 0.1 and a uniform electron density rig LR . In the 
following, primed quantities are measured in the reference frame comoving with the relativistic 
plasma blob, while unprimed quantities are generally referred to the stationary accretion-disk 
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frame. For simplicity, we assume that the relativistic electron distribution does not evolve 
during the outward-motion of the blob. This assumption leads to an overprediction of the 
contribution from the jet synchrotron mirror for which the accumulation of radiation along the jet 
trajectory is important. The synchrotron and Compton processes are calculated using 5- function 
approximations for the emission from relativistic electrons (see, e. g., Dermer et al. 1997). 

The important invariant with which to compute the energy density in accretion-disk radiation 
scattered by the BLR (denoted by the superscript 'ECC for EC radiation from Clouds) relates the 
angle-dependent photon number density between the comoving and the stationary frames, given 
by n' ph (e f , fi') = [r(l + /?r//)] -2 nph(e, fJ), where e = hu/(m e c 2 ) is the dimensionless photon energy. 
Following the treatment outlined by Bottcher & Dermer (1995), we find that the comoving- frame 
photon number density at height z 2 is given by 



/ECC, , x n* LR a T 2 f 'f NpjeuSli) 

n P h (^ z 2) = — y c — / dfl2 J dX2 — ^ — ' W 

..min ~,min 

where ot is the Thomson cross section, /j, 2 = cosa2 (cf. Fig. 1), e± = e' 2 /D2, and D 2 = T(l+(3rfi>2)- 
The term Nd(^i,^i) is the angle-dependent photon production rate of the disk in the direction 
^i = (Oi,<j>i)- In the derivation, note that the photon travels opposite to the direction for which 
a 2 is defined. The integration limits ^™ m / max anc j j™/™ i n equation (1) follow from the 
geometrical constraint that r- in < r < r out , and from the requirement x s > R' B , where x s = x 2 D 2 
is the distance of the location of the scattering event in the BLR to the center of the blob at the 
time of scattering in the BLR. The latter constraint represents the condition that there is no cold 
BLR material within the relativistic blob. 

The synchrotron reflection mechanism (GM96), where the jet synchrotron radiation is 
scattered within the BLR and reenters the jet at a later time and at a different location z 2 , is 
more complicated because light travel-time effects have to be carefully taken into account. The 
comoving-frame RSy photon number density at height z 2 is given by the expression 



n P h 



^raax ^max 

(e 2 ;^) = — y c — J dfi2 J dx z D i — ' ( > 



where e[ = e 2 /(DiD 2 ), D\ = [T(l — /?rMi)] _1 ) an( ^ Mi = cosai. N' sy {^i, z{) is the 
angle-dependent synchrotron photon production rate in the comoving frame at height z\, where 
Qi = (a[, </)[). In the derivation giving equation (2), we use the transformation relating the 
differential photon production rate between the stationary receiving and comoving frames, given 
by JVph(ei,fii) = DfN^n'J. 

The time of photon emission at location z\ is given by the causality condition 
z\ = z 2 — /3r(xi + x 2 ). Therefore, in addition to the geometrical constraints described in 
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the previous paragraph, the condition z\ > has to be met (note that the light travel-time 
constraint has been ignored in GM96). In particular, causality implies that the photon density 
due to the RSy process is 

n 'ph y = ) provided z 2 < z cr = r in ~ r in ~ ■ ( 3 ) 

For z CT < Z2 < n n , only radiation scattered from within a very small angle 
a max ^ p— i (2 2 //?r Hn) — 1 contributes to the rescattered photon density in the blob.[] 
Furthermore, since photons emitted at z = are at most a distance Az = ct{\ — /3p) ~ z/{2T 2 ) 
ahead of the blob when it reaches z, only a fraction Az/ArBLR of the radial Thomson depth of an 
extended cloud can scatter the synchrotron radiation from the blob. As viewed in the comoving 
frame, this means that unless r out 3> Ti?^, no extended photon halo (which contributes to the 
reflected radiation) will surround the blob as it passes through the BLR. 

A relativistically moving blob will receive a considerable flux of scattered synchrotron 
radiation only when it is passing through the BLR, i.e. for r; n <^ z <J r out . Assuming that 
ArBLR S> Ti?^, we find that the rescattered synchrotron photon energy density can be estimated 
as 

^4rV BL H^_(i-^)< r (4) 

Two powers of T come from Thomson scattering, and a third from length contraction and 
consequent density enhancement of scattering material as viewed in the comoving frame. 

In addition to the comoving-frame photon energy densities, u' ph = m e c 2 J °° de' e' n' ph (e'), we 
compute the total radiative power resulting from Compton scattering of the respective photon 
fields by calculating the total energy-loss rate of electrons due to Compton scattering. The decline 
of the Compton scattering cross section in the Klein-Nishina regime, e'7 ;> 1, is roughly accounted 
for by neglecting the energy-loss of electrons with Lorentz factor 7 which scatter photons of energy 
e' > I/7. The Compton power is 

! 72 1/7 

Kad = Vb-^cct t m e c 2 J dj 7 2 n e (7) J de' e' n ph (e'), (5) 
71 

where V' B denotes the blob volume. It is important to note that the different Doppler-factor 
dependences of the observed fluxes from different processes (especially the SSC and EC scattering 
processes; see Dermer 1995) are properly taken into account. 



1 The small solid angle illuminated by the RSy mechanism in comparison with the large solid angle subtended by 
the BLR clouds in the ECC process could explain why no varying Lya emission correlated with the highly variable 
UV continuum was detected from 3C 279 (Koratkar et al. 1998). 
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3. Numerical Results and Discussion 

Photon energy densities and radiative powers due to different processes were numerically 
calculated as a function of the distance of the blob from the accretion disk for a wide range 
of parameters. Figure 2 shows an example of our series of simulations, using parameters 
representative of a flat-spectrum radio quasar (FSRQ) and a BL Lac object (BL). Here we let 
p = 3 (FSRQ; solid curves) and p = 2.7 (BL; shaded curves). This choice yields fairly flat vF v 
spectra which can be compared with the peaks of the broadband vF v spectral energy distributions 
of blazars. The electron density and blob radius are chosen so that the resulting total luminosities 
are in accord with typical values of the apparent luminosities of FSRQs and BLs, and so that 
< 100 GeV 7 rays are not absorbed by 7-7 pair production on the synchrotron photons intrinsic 
to the source. 

For the generic FSRQ, we assume that the BLR has a radial Thomson depth tt,blr = 0-2 and 
occupies a spherical shell located between 0.05 and 0.5 pc from the central engine. For simplicity, 
we assume that the accretion disk radiates isotropically with a luminosity of 10 46 erg s _1 . This 
choice of parameters gives a Compton power which is ~ 10 times greater than the synchrotron 
power. Figure 2 shows that although the RSy radiation energy density u' RSy is larger than u' Sy 
when the blob is located within the BLR, the RSy radiative power is about equal to the SSC 
power. This is because the backscattered jet synchrotron photons are boosted in energy by a factor 
of r 2 relative to the energy of the synchrotron photons in the comoving frame; thus Klein-Nishina 
effects reduce the radiative power resulting from Compton scattering of RSy radiation more 
strongly than for the SSC radiation. 

In our simulations using plausible parameters compatible with the assumed spherical shell 
geometry of the BLR, we find that the Compton power due to the synchrotron mirror mechanism 
is at most comparable to the SSC power. The efficiency of this process is improved when T S> 10 
and a very narrow (Afblr <S Hn/(2r 2 )), Thomson-thick BLR is located very far (r- m 3> r 2 i?' B ) 
from the central engine. The synchrotron reflection process might therefore operate in BLR 
clouds which are thought to surround Seyfert AGNs and FSRQs. BLR clouds, as understood 
through photoionization models (see, e.g., Wandel 1997 and references therein), consist of dense 
(n ~ 10 10 — 10 11 cm), Thomson-thick (r ~ 1 — 10) regions covering a small (~ 10%) fraction of 
the central engine. For such a model to be feasible, however, the conditions regarding duty cycle 
and power outlined by Dermer & Chiang (1998) must be met. We defer presentation of results in 
the regime R' B ;> AreLR/r to future work; the shell geometry used here is excessively artificial in 
this limit. 

The absence of strong emission lines in X-ray selected BLs and (to a lesser extent) 
radio-selected BLs suggests that the BLR is considerably more dilute in BLs than in FSRQs, and 
that BLs have mean Thomson thicknesses tt,blr "C 0.1 (see Scarpa & Falomo 1997 and references 
therein). (On the other hand, the strength of the central ionizing photon source might be much 
less in BLs than FSRQs.) Superluminal motion observations also indicate that typical values of T 
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for BLs lie in the range between ~ 3 and 7 (see the review by Urry & Padovani 1995). The ability 
of the synchrotron reflection process to produce gamma-ray flares therefore seems more difficult 
in BLs than in FSRQs, yet TeV flares often coincident with X-ray flares have been detected from 
three BL Lac objects (e.g., Punch et al. 1992, Macomb et al. 1995; Catanese et al. 1997, 1998). 
For the assumed BL parameters in Figure 2, the synchrotron reflection flare could hardly be 
detected. 

If an accretion disk steadily radiates photons (see Bottcher & Dermer 1995 for a treatment of 
time- variable disk radiation), then light-travel time effects can be neglected. The accretion disk 
supplies an abundant supply of soft photons, which can enter the jet directly (ECD) and after 
being scattered by the BLR (ECC). The comoving photon energy densities from the ECD process 
can dominate that from the ECC and synchrotron processes when z <J 10~ 2 pc, but declines oc z^ 
and oc z~ 2 farther out (see Dermer & Schlickeiser 1993; Bottcher, Mause & Schlickeiser 1997). 
The ECC photon energy density increases slowly with z when z < r- m , and begins to decrease for 
j"i n <; z <; r out . Outside the BLR, when z ;> r out , the energy density asymptotically approaches the 
limiting behavior u' ECC oc z~ 2 . 

When z <^ r m , the ECC process dominates over the SSC process provided that 



7"T,BLR 2 cB z T B ( P~ 1 . 7 2 ~7l 

rf n * 3 Ld r 2 { 3-p> 7l 1 - p -7 2 1 - J " W 
where tb = R'b n c,jct o"T is the Thomson depth of the blob. Using the parameters adopted 
in Figure 2 but letting tt,blr> r im and L vary, we find that the ECC photon energy 
density dominates the synchrotron photon energy density in the comoving blob frame when 

1/2 1/2 1/2 1/2 

r in(pc) ^ 0-4T T ' BLR L 4 g anc l rj n (pc) ^ 0.1 Tt BLR L 4 4 f° r the FSRQ and BL parameters, 
respectively, where L n = Ldi s k/(10 n ergs -1 ). It should be noted that Klein-Nishina effects are not 
included in estimate (6). 

The bottom panel of Figure 2 illustrates the time profiles of a flare calculated using the FSRQ 
parameters for the ECC (thick solid curve) and RSy (thick dot-dashed curve) processes. Note that 
the observer's time element is linearly related to z for a blob moving with constant velocity. The 
ECC process gives a fast-rise, power-law-decay-type light curve, and the RSy mechanism gives a 
gradual rise of the 7-ray flux and a sharp drop as the blob leaves the BLR. A flare produced by 
the RSy process could be identified by a rapid decline of 7-rays which is not accompanied by a 
corresponding decrease of the synchrotron emission. For comparison, we also sketch a flare time 
profile produced by the ECD process, and time profiles produced by sweeping energization of the 
blob. In this process, the bulk kinetic energy of the outflowing plasmoid is converted into internal 
nonthermal particle energy by sweeping up BLR material (see, e.g., Panaitescu & Meszaros 1998; 
Dermer & Chiang 1998; Chiang & Dermer 1998). The pair of light solid and dot-dashed curves 
illustrate ECC and RSy flares, respectively, are modeled assuming that the nonthermal electron 
energy is proportional to the amount of swept-up matter which is then added to a nonthermal 
lepton distribution accelerated at the base of the jet. Blob deceleration is assumed to be negligible 
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here. The upper curves of the two pairs are modeled assuming no radiative losses, and the lower 
curves of the two pairs illustrate the effects of radiative losses on the nonthermal leptons by 
crudely multiplying the upper curves by a decaying exponential with a 1.5 x 10 18 cm decay length. 

The model RSy light curves, with and without sweeping energization, are similar to 7-ray 
light curves observed in the 1991 and 1996 flares of 3C 279 (Hartman et al. 1996; Wehrle et 
al. 1998). This may be a consequence, however, of the highly idealized BLR geometry used in 
the calculation. More symmetical flaring profiles observed from PKS 0528+134 (Collmar et al. 
1997), PKS 1622-297 (Mattox et al. 1997), and PKS 1406-076 (Wagner et al. 1995) might be 
more easily explained by the ECD or ECC processes. The declines of the X-ray and optical fluxes 
correlated with the EGRET 7-ray fluxes in the February 1996 3C 279 and the 1406-076 flares 
could, however, rule out the RSy mechanism since the synchrotron component is not directly 
affected by the reflection process. The inclusion of electron energy evolution and the relaxation of 
the assumption of a constant velocity blob must be treated to strengthen such conclusions. Future 
flare modeling must treat the passage of the jet through the BLR clouds, which themselves are in 
Keplerian motion around the central black hole. The passage of a jet through such a region will 
display a complicated signature when monitored by different telescopes with different sensitivities 
and imaging capabilities, which can only be decoded when full account is taken of the processes 
considered here. The efficiency of the different Compton-scattering scenarios, including the RSy 
mechanism, in this more realistic system is presently under investigation by the authors. 

4. Summary 

We have analyzed different mechanisms for producing 7-ray flares which involve relativistic 
leptons in blazar jets which Compton-scatter soft photons. We presented a detailed study of 
the jet synchrotron reflection process (Ghisellini & Madau 1996), where the beamed synchrotron 
radiation emitted by a relativistic plasma blob is scattered back into the jet by BLR material, and 
determined conditions under which this mechanism can operate efficiently compared to competing 
mechanisms. Our calculations show that the efficiency of the synchrotron reflection mechanism is 
improved for very narrow (AreLR ^ nn/2r 2 ) Thomson-thick tblr ^ 0.1 scattering regions located 
far (r- m ^0.1 pc) from the central accretion disk. We argued that the synchrotron reflection model 
may have serious difficulties explaining 7-ray flares from BL-Lac objects. Model flare time profiles 
for various Compton-scattering processes, which in addition include effects of energization of the 
blob by sweeping up circumnuclear material, were illustrated. The slow-rise, rapid-decay time 
profile seen in several EGRET flares from FSRQs could be realized by the RSy process, though 
an accompanying flare from the ECC and ECD processes would also be seen unless the blob 
only becomes energized after it enters the BLR. A more detailed study of flare production which 
includes the reflected synchrotron mechanism and sweeping energization in a system where the 
BLR consists of small, Thomson thick clouds is an important subject for further study. 



-8- 



MB acknowledges support by the German Academic Exchange Service (DAAD) and by 
NASA grant NAG 5-4055. The work of CD is supported by the Office of Naval Research and the 
Compton Gamma Ray Observatory Guest Investigator Program. 

REFERENCES 

Blandford, R. D., & Konigl, A. 1979, ApJ, 232, 34 
Blandford, R. D., & Levinson, A. 1995, ApJ, 441, 79 
Bloom, S. D., & Marscher, A. P., 1996, ApJ, 461, 657 

Bloom, S. D., Thompson, D. J., Hartman, R. C, & von Montigny, C, 1997, in 4th Compton 
Symposium, ed. C. D. Dermer, M. S. Strickman, & J. D. Kurfess (New York: AIP), 1262 

BSttcher, M., & Dermer, C. D., 1995, A&A 302, 37 

Bottcher, M., Mause, H., & Schlickeiser, R., A&A, 324, 395 

Catanese, M., et al. 1997, ApJ, 487, L143 

Catanese, M., et al. 1998, ApJ, in press 

Chiang, J., & Dermer, C. D. 1998, ApJ, submitted 

Collmar, W., et al., 1997, in Proc. of the XXV. ICRC, Vol. 3, 101, ed. M. S. Potgier, B. C. 
Raubenheimer, & D. J. van der Walt 

Dermer, C. D., 1995, ApJ, 446, L63 

Dermer, C. D., & Chiang, J., 1998, New Astronomy, 3, 157 
Dermer, C. D., & Schlickeiser, R., 1993, ApJ, 416, 458 
Dermer, C. D., Sturner, S. J., & Schlickeiser, R., ApJS, 109, 103 
Ghisellini, G., & Madau, P., 1996, MNRAS 280, 67 
Hartman, R. C, et al. 1996, ApJ, 461, 698 

Hartman, R. C, Collmar, W., von Montigny, C, & Dermer, C. D., in The Fourth Compton 
Symposium, ed. C. D. Dermer, M. S. Strickman, & J. D. Kurfess (New York: AIP), 307 

Hummel, C. A., Krichbaum, T. P., Witzel, A., et al., 1997, A&A, 324, 857 

Koratkar, A., Pian, E., Urry, C. M., & Pesce, J. E. 1998, ApJ, 492, 173 

Macomb, D. J. et al. 1995, ApJ, 449, L99; (e) 1996, 459, Llll 

Maraschi, L., Ghisellini, G., k Celotti, A., 1992, ApJ, 397, L5 

Marscher, A. P., & Gear, W. K., 1985, ApJ, 298, 114 

Marscher, A. P. 1990, in Parsec-Scale Radio Jets, ed. J. A. Zensus & T. J. Pearson (New York: 
Cambridge University Press), 236 



-9- 



Mattox, J. R., et al., 1997, ApJ, 481, 95 

Mukherjee, R., Bertsch, D. L., Bloom, S. D., et al., 1997, ApJ 490, 116 

Panaitescu, A., & Meszaros, P. 1998, ApJ, 492, 683 

Protheroe, R. J., & Biermann, P. L. 1997, Astroparticle Physics, 6, 293 

Punch, M., et al., 1992, Nature 358, 477 

Scarpa, R., k Falomo, R. 1997, A&A, 325, 109 

Schlickeiser, R., 1996, Space Sci. Rev., 75, 299 

Shakura, N. I., Sunyaev, R. A., 1973, A&A, 24, 337 

Sikora, M., Begelman, M. C, & Rees, M. J., 1994, ApJ, 421, 153 

Steffen, W., et al., A&A, 302, 335 

Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 

Wagner, S. J., et al. 1995, ApJ, 454, L97 

Wandel, A. 1997, ApJ, 490, L131 

Wehrle, A., et al. 1998, ApJ, in press 



This preprint was prepared with the AAS IATjTjX macros v4.0. 



-10- 




Fig. 1. — Geometry used for the computation of jet radiation rescattered into the jet trajectory. 
Radiation is emitted at distance z\ from the accretion disk at z = 0, is rescattered in the BLR at 
a distance r from the central engine, and reenters the jet at location Z2- 
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z [cm] 

Fig. 2. — Energy densities (upper panel), radiative powers (middle power), and model light cuves 
(bottom panel) for synchrotron and Compton processes, plotted as a function of distance z of the 
blob from the central engine. 'ECC denotes accretion-disk photons scattered from the BLR into the 
blob, 'RSy' denotes reflected synchrotron photons, 'ECD' denotes accretion-disk photons entering 
the blob directly, and '+E' curves illustrate the effects of sweeping energization on the blob as it 
passes through the BLR. Standard parameters used for a typical FSRQ (solid curves) are T = 10, 
7l = 10 3 , 72 = 2 • 10 5 , p = 3.0, n ei j et = 20 cm" 3 , B = 0.91 G, R' B = 2 • 10 16 cm, r in = 0.05 pc, 
^out = 0.5 pc, tt,blr = 0-2, and L^isk = 10 46 erg s _1 . Standard parameters used for a typical BL 
(shaded curves) are the same as for the FSRQ case except that 72 = 10 6 , p = 2.7, R' B = 10 16 cm, 
B = 0.99 G, tt,blr = 0.02, and Ldisk = 10 44 erg s _1 . Bottom panel uses a linear scale and FSRQ 
parameters. 



